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Abstract. 

Turbulent combustion is three-dimensional. Turbulence in a Type la su- 
pernova is driven on large scales by the buoyancy of burning products. The 
turbulent cascade penetrates down to very small scales, and makes the rate 
of deflagration independent of the microphysics. The competition between 
the turbulent cascade and the freeze-out of turbulent motions due to stellar 
expansion determines the largest scale participating in the cascade. This 
sets the bulk rate of a deflagration in a supernova. The freeze-out limits the 
bulk rate of deflagration to a value that makes a powerful explosion impos- 
sible. Two-dimensional simulations cannot capture these essential elements 
of turbulent combustion, even in principle. 

A powerful delayed detonation explosion can take place if the burning 
makes a transition to a detonation. A deflagration to detonation transition 
(DDT) can occur in a layer of mixed cold fuel and hot burning products 
created either inside an active turbulent burning region by a high intensity 
of turbulence, or by mixing cold fuel with ashes of a dead deflagration front 
during the global pulsation of a star, or by both. 



1. Introduction 



We discuss three major topics: 

1. Turbulent combustion in a uniform gravitational Held. We summarize the 
results of three-dimensional numerical modeling of buoyancy-driven turbu- 
lent burning in supernovae and in terrestrial environments and point out 
the possibility of terrestrial experiments designed specifically for modeling 
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terrestrial turbulent combustion and understanding combustion in super- 
novae. 

2. Turbulent combustion in an expanding star. We describe the Grst three- 
dimensional fluid dynamics simulations of the deflagration of an entire 
carbon-oxygen white dwarf exploding as a Type la supernova. The results 
clearly point to the important role of the expansion freeze-out of turbulent 
motions. 

3. Deflagration to detonation transition (DDT). Based on a theory of de- 
flagration to detonation transition (DDT) in an unconfined environment, 
and on the estimates of the turbulent velocity obtained from the three- 
dimensional simulations, we analyze the possibility of the DDT in super- 
novae, and its relevance to the (pulsating) delayed detonation models of 
Type la supernova events. 

Limitations of the current state of numerical modeling, and important prob- 
lems for the future are discussed in the conclusion. 

Nomenclature 



R, 


- stellar radius 




- expansion velocity 


9 


- gravitational acceleration 




- flame radius 


Lf 


- freeze-out length 




- expansion time-scale 


L 


- column width 


Si 


- laminar flame speed 


5i 


- laminar flame thickness 


F 


- Froude number 


U' 


- turbulent velocity 


Ag 


- Gibson scale 


Ac 


critical wavelength for 


a > 1 


- density ratio across 




the onset of turbulent burning 




laminar flame front 


D 


- detonation wave speed 


Sd 


- detonation wave thickness 




- detonation cell size 


Afe 


- detonation kernel size 


Li 


- detonation initiation length 


Ui 


- detonation initiation velocity 



2. Turbulent combustion in a uniform gravitational field 

The issue of turbulent combustion in a uniform gravitational field singles 
out one important element of supernovae deflagration - propagation of the 
flame against gravity. The turbulence is not imposed externally, as in many 
combustion experiments and theoretical studies, but is driven internally 
by the buoyancy of burning products, as in supernovae. This problem can 
be studied both by direct numerical simulations and experimentally. Using 
the results of studies of this restricted problem, various subgrid numerical 
models for turbulent combustion can be calibrated and then applied to 
supernovae. 



3 



Consider a square column of size L filled with a premixed gaseous fuel 
in a uniform gravitational field g. The fuel is initially in the state of hydro- 
static equilibrium. After the fuel is ignited at the bottom, a flame begins to 
propagate vertically. Because the propagation is against gravity, the flame 
becomes turbulent due to the RT instability ^. 

For uniform conditions, the flame would eventually reach a steady-state 
burning velocity, St, which is statistically independent of time. The general 
functional dependence of St on the parameters of the problem follows from 
dimensional arguments: St = sfgL f{a, F), where F = Sf/gL is the Froude 
number which characterizes the relative importance of laminar burning at 
speed Si and the R-T instability to the fluid flow on the largest scales. 
When gravity dominates (F 0), the turbulent flame velocity becomes 
independent of the details of small scale motions, and of the microphysics 
(nuclear reaction rates, etc.) of burning. The turbulent flame speed becomes 
a function of g, L, and a only, 

St = VgLfia) c c [gL—^) . (1) 

Numerical simulations confirm this scaling law and give C c± 0.5 for a 
square column with periodic boundary conditions [5] [6] . 

The mechanism underpinning the scaling law is the competition of two 
effects: surface creation by stretch, and surface destruction due to self- 
intersections of the highly distorted flame. In a steady state, these effects 
balance each other, and 5"^ becomes independent of Si. The flame regulates 
itself so as to create more surface if Si is small, and creates less surface 
if Si is large [5]. The scaling law suggests a linear relation between the 
turbTilcnt flame speed and the velocity of turbulent motions on the largest 
scale, St ^ U'{L). It is clear that the self-regulating mechanism operates in 
all situations when the turbulence is driven on large scales, and cascades 
down to small scales. 

For this self-regulation to take place, the turbulent scale where the tur- 
bulent velocity equals the laminar flame speed, also called the Gibson scale 

laminar flame front propagating against gravity is unstable with respect to pertur- 
bations of all wavelengths exceeding ~ 10 — 100<5;, where 5i is the laminar flame thickness 
[1][2]. It is nonlinearly stabilized due to the difference in the behavior of concave and 
convex parts of the flame surface and formation of cusps [3] for a wavelength of pertur- 
bations A < Ac — — — ^^ryi where 5"; is the laminar flame speed, g is the acceleration 
of gravity, and q > 1 is the ratio of densities ahead and behind the flame front [5]. For 
supernovae conditions, Ac ~ 10"' — 10^ cm ^ 5i, depending on density and position of 
the flame in the star. On larger scales, the flame is truly turbulent. The transition to 
turbulence is facilitated by the development of secondary instabilities, and occurs on a 
time scale of ~ 2 — 3 RT characteristic times. 
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Figure 1, Regimes of combustion on the f'^/F — L/5i plane. Line (1) - f'^ /F = 25. Line 
(2) \g = 2(5;. Laminar and weakly wrinkled flanielet regimes are below the line (1). 
The multiple flame sheet or corrugated flamelet regime is between the lines (1) and (2). 
The scaling analysis of §2 applies to this regime. The distributed flame regime is above 
line (2). The possibility of initiating a detonation inside an active turbulent flame front 
(§3) exists in this regime. Circle - conditions of the terrestrial experiment discussed in 
§2. Triangle - conditions for turbulent burning in the middle of a supernova with p = 10* 
g/cm^ , at the scale L = 2 x 10® cm. 

[7], can be expressed as: 

(2, 

must be larger than the laminar flame thickness 6i. This self-regulated 
regime can be classified as a corrugated fiame regime [7] , or as a regime of 
multiple flame sheets [8]. 

One example of a terrestrial system well suited for experimental study 
of buoyancy-driven turbulent burning is a lean methane flame diluted to an 
equivalence ratio (j) ~ 0.55 — 0.6. The laminar flame velocity is Si ~ 10 cm 
s~^, and the flame thickness is (5; ~ 1 mm [14]. The regime of the scaling 
law of equation (1) can be achieved in a vertical tube with L ~ 30 — 50 cm, 
and height of several meters [6]. This regime is shown in Figure 1 which 
shows the plane spanned by the parameters /^/F and L/5i. The circle shows 
the condition predicted for the methane experiments and the triangle shows 
burning on the scale L = 10^ cm in the middle of a supernova {g = 1.9 x 10^ 
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cm/s, Si = 10^ cm/s). Experiments to check the predictions of the scahng 
law are in progress. 

3. Turbulent combustion in an expanding star 

The uniform gravitational field approximation could be valid on spatial 
scales less then the mean radius of the flame Rf, and on time scales less 
than the time scale of the expansion Te ^ dt/dln R^, where i?* is the stellar 
radius. The main effect of expansion is to freeze turbulence on the scales 
where the turbulent velocity generated by the RT instability is comparable 
to or less than the expansion velocity on those scales. We expect this freeze- 
out to occur on a scale Lf where the turbulent burning velocity given by 
equation (1) becomes comparable with the expansion velocity on this scale 
[5], i.e, where 

St{Lf)^Ue{Lf)=Lf/Te . (3) 

To get an estimate of we use for the expansion rate the results of a 
one-dimensional simulation with the average turbulent flame speed given 
by equation (1) with L equal to the flame radius (assuming no turbulence 
freeze-out or effects of spherical geometry) [5]. This gives Lj ~ (a few) x lO'^ 
cm, less than the mean radius of the flame Rf which exceeds (a few) x 10^ 
cm during the final stages of deflagration, and clearly shows that freeze-out 
is important. Using equation (1), we then estimate the effective deflagration 
speed as 




where is the sound speed in the stellar matter. In conditions typical of 
an exploding white dwarf, we do not expect the turbulent burning speed 
to exceed a few per cent of the sound speed. This is not enough to cause a 
powerful explosion even in the idealized one-dimensional picture where the 
front is spherical, and all matter passing through the front is burned. 

An additional effect that further limits the rate of deflagration is a 
deviation from the steady-state turbulent burning regime. A certain time, 
2± \/St{X), is required for a turbulent flame to reach a steady-state on the 
scale A. This time is larger for larger scales. Scales of the order oi Rf might 
never reach a steady-state. This can be also expressed in a different way. 
The thickness of a steady state turbulent front in a vertical column of size 
L is 2± 2 — 3L [5] [6]; but for L 2± there is not enough space in a star to 
fit a steady-state turbulent flame inside the star. 

Some results of a three-dimensional numerical simulation of the entire 
white dwarf exploding as a supernova are presented in Figures 2 and 3 (see 
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Figure 2. Results of a three-dimensional simulation of a turbulent deflagration explosion 
of a IAMq carbon-oxygen white dwarf (§3). The central density (dashed line) and the 
amount of burned matter (in per cent) as a function of time elapsed since ignition in the 
center are shown. 



[5] for details). In these simulations, equation (1) has been used for the 
turbulent flame velocity on scales not resolved numerically. For regions of 
the "average" flame front not oriented "upwards" against gravity this for- 
mula most probably overestimates the local turbulent flame speed. Despite 
this, only 2± 5% of the mass has been burned by the time the model has 
expanded and quenched the flame, and the white dwarf has not become 
unbound. These results show that burning on large scales does not reach a 
steady state. Big blobs of burned gas rise and penetrate low density outer 
layers, whereas unburned matter flows down and reaches the stellar center. 
The flame was not able to close the gaps between the blobs and make the 
flame surface quasi-spherical. Even if this had occurred, however, the speed 
of the deflagration would have been exactly as described by equation (1) 
with L given by equation (3), and this is not enough to unbind the model. 

The model experienced an almost complete overturn. This has obviously 
important implications for nucleosynthesis and may cause an element strat- 
ification incompatible with observations if composition inhomogeneities are 
not smeared out during the subsequent detonation stage of burning. 

The results of 3D modeling indicate that the deflagration alone is not 
sufficient to cause an explosion. To make a powerful explosion, the deflagra- 
tion must somehow make a transition to a detonation (delayed detonation 
model). This idea is also supported by observations of SN la, events that 
are best described in the framework of the delayed and pulsating delayed 
detonation scenarios ([18], [19], [20]). 
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Figure 3. The distribution of burned inside a 1.4A/0 white dwarf resulting from ~ 1.5 
s of turbulent combustion for the simulation shown in Figure 2. 

4. Deflagration to detonation transition (DDT) 

In terrestrial conditions, DDT is often caused by interaction of a turbu- 
lent flame with obstacles or boundaries. DDT in unconfined conditions is 
possible, but more difficult to achieve. Here we apply the theory of DDT 
in terrestrial unconfined conditions proposed in [21] to a DDT in super- 
novae. The theory is based on two major assumptions: {i) The detonation 
arises from a non-uniform explosion of a region of a fuel with a gradient 
of induction time via the Zeldovich gradient mechanism [22] [23]. This re- 
gion is created by mixing of fresh fuel and hot products of burning, and 
is characterized by the oppositely directed gradients of temperature and 
composition. There exists a minimum size of the region capable of generat- 
ing a detonation, Lj. In [21] we argue on the basis of numerical simulations 
and terrestrial experiments that Li can be correlated with the experimental 
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values of the detonation kernel size Ik (the minimum length to propagate a 
steady-state detonation) as ~ (1 — 4)/^ or with the detonation cell size 
Ic (the size of a single detonation "cell" defined by the loci of intersecting 
triple points in the wake of a detonation front) as Lj ~ 10 — 40/c, or with 
the theoretical value of the one-dimensional detonation wave thickness 5d, 
as Lj ~ W^Si- The exact coefficients in these relations depend on the de- 
tails of reaction kinetics and equation of state, (ii) A level of turbulence Ui 
required to mix the products and fuel inside a region of size Lj follows from 
the requirement that the Gibson scale given by equation (2) be comparable 
to the laminar flame thickness Si, or 

Ui = KSi {Li/5ifl^ (5) 

where the coefficient K ~ 1 — 10 describes the ability of the flame to stabilize 
itself against turbulent motions on scales ~ (5/. At this level of turbulence, 
a sharp interface - a laminar flame front - between the fuel and products 
cannot survive turbulent stretch. The flame front can be locally broken 
and extinguished, and the mixing of fuel and products can occur on scales 
» Si. 

The one-dimensional thickness of a detonation front in a carbon-oxygen 
fuel (C+C reaction zone) ranges from ~ 10~^ cm at densities ~ 10^ g/cm^ 
to ~ 10 cm at ~ 10''^ g/cvc? [24]. From this we can roughly estimate Lj ~ 
10 — 10^ cm for the same density range. The two-dimensional structure 
of a detonation front has been computed for the density 3 x 10^ g/cvc? 
in [29]. The detonation cell size was found to be ~ 10 cm, which gives a 
slightly smaller value of Li ~ 10^ cm than the one-dimensional estimate 
Li ~ 10^ cm. It is, however, clear that the estimates are correct within 
an order of magnitude. We also notice the ratio of thicknesses of the one- 
dimensional detonation and a laminar flame in a degenerate CO mixture is 
approximately Si/Sd ^ L)/Si ~ 10~^. 

With these data in mind, we can rewrite equation (5) in numbers as 

El ^ 10-2 ( 5 X 10^cm/s \ /^^X 

The laminar flame speed ranges from 10^ to 10'* cm/s for densities 
10^ — lO'^ g/cm^, respectively [26]. Comparing equations (4) and (6) we see 
that it might be possible to trigger a detonation inside an active turbulent 
flame front when burning reaches densities ~ 10^ g/cm^. 

If an active turbulent deflagration fails to make the transition to a 
detonation, stellar expansion will eventually quench nuclear reactions at 
densities less than ~ 10^ g/cm^. The star will experience a pulsation and 
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collapse back. During the expansion and contraction phases of the pulsa- 
tion, the high entropy ashes of the dead deflagration front will mix with 
the the fresh low entropy fuel again to form a mixture with temperature 
and composition gradients. Mixing will be facilitated during the contraction 
phase because of the increase of turbulent motions due to the conservation 
of angular momentum. An estimate of the mixing region formed during 
pulsation is ~ 10^ — lO''' cm, much larger than Lj. As soon as this mixture 
returns to high enough densities ~ 10^ and re-ignites, the detonation will 
be triggered ([19], [25]). 

The theoretical conclusion is that DDT should occur in a carbon-oxygen 
white dwarf one way (delayed detonation) or another (pulsation delayed 
detonation) at densities ~ 10'' g/cm^, but at present it is difficult to say 
which mode is realized more frequently. The estimate of the transition den- 
sity at which detonation must be triggered is consistent with the estimate of 
the transition density derived from fitting one-dimensional delayed detona- 
tion and pulsation delayed detonation models with observations of SN la. A 
recent success of the delayed detonation model is the explanation of the cor- 
related photometric and spectroscopic diversity of SNIa (([20], [15], [16], [27], 
[28], [17]; see also contributions by Hoflich et al. in these proceedings). 

5. Conclusions 

Despite substantial progress in hydrodynamical modeling and understand- 
ing of SNIa, there are many important problems yet to address. A few 
important ones are: 

To determine the exact deflagration velocity inside a carbon-oxygen white 
dwarf as a function of initial composition and initial conditions. This re- 
quires high resolution 3D modeling in order to correctly represent the com- 
petition between the turbulent cascade and turbulence freeze-out. 

To understand the mode(s) of deflagration to detonation transition in a 
white dwarf. This would require both large scale modeling of the pulsation, 
and small scale modeling of mixing processes in active and dead deflagration 
fronts. 

Both efforts will require new numerical techniques to provide a resolution 
of many orders of spatial scales in 3D. One such a technique, an Eulerian 
Tree adaptive mesh refinement, has been recently developed in Texas and 
successfully applied to the 3D problem of stellar disruption ([30]). It is 
being applied to both terrestrial and astrophysical combustion problems 
(in progress). 
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